A Supernova Factory in the Merger System Arp 299 



Susan G. Neff 

NASA 's Goddard Space Flight Center, Code 681, Greenhelt, MD 20111 

nef f Ostars .gsf c .nasa.gov 

^". 

^ ■ James S. Ulvestad 

National Radio Astronomy Observatory, P. 0. Box O, Socorro, NM 81801 

julvestaOnrao . edu 



o 

(N 



1—5 

oo 



o 



Ph' 



Stacy H. Teng 
Department of Astronomy, University of Maryland, College Park, MD 20142 

stacytOastro . umd . edu 



(N 

^^ , We have imaged the nearby galaxy merger Arp 299 at arcsecond and milharcsecond resolution, 



ABSTRACT 



using both the Very Large Array and the Very Long Baseline Array. The large-scale radio 

emission from the merger contains 5 bright, compact radio sources embedded in diffuse emission, 

1-^ ' with diameters less than 200 pc. Supernova rates of 0.1 to 1 per year are required to produce 

the VLA-detected radio emission in these sources. Two of the compact VLA radio sources, 



O ' designated Source A and Source D, also have been detected and imaged at milharcsecond scales. 

^ , Source A, which is associated with the nucleus of one of the merging galaxies, contains five 

^ ' milliarcsecond-scale sources, each with a radio power between 100 and 1000 times that of the 

Galactic supernova remnant Cassiopeia A. Four of these have flat or inverted spectra and appear 
to be young supernovae. Three of the VLBI-scale sources are located within 10 pc (projected) 
^^ , of one another, and two are separated by less than 3 pc, indicating that they all may be within 

Vh ' the same super starcluster or complex of such clusters. The brightest VLBI-scale source, AO, 

^ has an extremely inverted spectrum, with a > +2 at gigahertz frequencies. It seems to be the 

youngest supernova, which has not yet broken out of its circumstellar shell. The milharcsecond 
radio sources within Source A appear to constitute a supernova factory, confirming the presence 
of an extreme starburst that peaked at least a few million years ago. 

Subject headings: galaxies: individual (Arp 299, NGC 3690) — galaxies: interacting — galaxies: starburst 
— stars: clusters 

1. Introduction 

Arp 299 (NGC 3690, Mrk 171)^ is the "origi- 
nal" starburst galaxy (Gehrz, Sramek, & Weed- 



^The Arp 299 merger is sometimes referred to as NGC 3690 
(west) and IC 694 (east), although the designation of vari- 
ous merger components and nearby companions in the lit- 
erature is inconsistent (see Yamaoka et al. (1998) for a dis- 
cussion); in this paper we refer to the whole system as 
Arp 299. 



man 1983) and an obvious merger system that 
has been studied extensively at many wave- 
lengths. The galaxy distance is 41 Mpc for 
Ho = 75 km s^^ Mpc-^ (TuUy 1988), where 1" 
corresponds to 200 pc. The merger in Arp 299 is 
similar in age to the "Antennae" ; both systems 
are considered to be in the "early" merger stage, 
slightly more than one crossing time (> 10^ yrs) 
before coalescence of the nuclei is expected (Keel 
& Wu 1995; Read & Ponman 1998). An active 
starburst in Arp 299 is indicated by the high fre- 
quency of recent optically discovered supernovae 
in the galaxy, including SN 1992BU (Van Buren et 
al. 1994), SN 1993G (Forti 1993), SN 1998T (Li, 
Li, & Wan 1998; Yamaoka et al. 1998), and SN 
1999D (Qiu, Qiao, & Hu 1999). Lai et al. (1999) 
used 2.2-/xm adaptive optics imaging of Arp 299 
and detected what they believe to be many "su- 
per star clusters" (SSCs) that might include very 
high mass stars. Hubble Space Telescope (HST) 
FOG and NICMOS images (Meurer et al. 1995; 
Alonso-Herrero et al. 2000) reveal a population of 
young stellar clusters at ultraviolet, optical and 
near-infrared wavelengths. 

At visible wavelengths, the disks of the interact- 
ing galaxies appear to overlap, but the individual 
nuclei still are separate (~ 4 kpc) and distinguish- 
able (See Figure 1). Unlike the "Antennae", which 
has two clearly distinguished tidal tails, Arp 299 
has a single tidal tail, indicative of a prograde- 
retrograde interaction in the past (^700 Myr ago). 
The most recently initiated interaction also ap- 
pears to be prograde-retrograde (Hibbard 1997), 
a spin geometry which drives most of the gas to- 
wards the galaxy center and maintains star for- 
mation at an elevated level for an extended time. 
Most of the molecular gas in the system is con- 
centrated near the nuclei, particularly the east- 
ern (retrograde) nucleus (Sargent & Scoville 1991), 
but CO is also found in other regions such as the 
bridge between the nuclei and the base of the tidal 
tails (Aalto et al. 1997; Casoh et al. 1999). 

Arp 299 's infrared luminosity {Lm > 3 x 
IO^^Lq) is nearly in the Ultraluminous IR Galaxy 
(ULIRG) category (Hibbard 1997). The eastern 
nucleus, hereafter source A (Gehrz et al. 1983), 
is the brightest component in the near-infrared, 
accounting for ~50% of the total infrared lumi- 
nosity of the system (Alonso-Herrero et al. 2000). 
There are numerous other point-like near-infrared 



sources in Arp 299, many of them in very dusty re- 
gions of the system (Alonso-Herrero et al. 2000). 
Many, but not all, of the infrared objects corre- 
spond to sources detected at radio wavelengths in 
VLA^ images (Gehrz et al. 1983; Zhao et al. 1997). 
Numerous H II regions populate the system near 
star-forming regions and it appears that star for- 
mation has been occurring at a high rate for about 
10 Myr (Alonso-Herrero et al. 2000). 

The new and archival radio observations pre- 
sented here are part of a larger project studying 
the star formation rate of mergers, to determine 
at which stage of the merging process the maxi- 
mum star formation occurs. This paper focuses on 
the properties of the five strongest radio sources in 
Arp 299, investigating them on scales of ~ 100 pc 
and of <1 pc. The population of weaker radio 
sources and the diffuse radio emission in the sys- 
tem will be the subject of a separate paper (Neff, 
Gampion, and Ulvestad, in prep). 

2. VLA Observations, Archival Data, and 
Imaging 

New observations of Arp 299 were made us- 
ing the VLA (Thompson et al. 1980) at 4.9 and 
8.4 GHz; these wavelengths were chosen because 
they are the most sensitive bands on the VLA. The 
new 4.9 and 8.4 GHz observations were made on 
2000 October 24, 2001 March 19, and 2002 April 
29. In 2000 October, the system was observed at 
4.9 GHz with the A configuration. Eight-minute 
integrations of the galaxy were alternated with 1- 
minute scans of the phase calibrator Jl 1284-5925 
over a period of four hours; two 2-minute observa- 
tions of 3G 286 were used for amplitude calibra- 
tion. Arp 299 was observed at 8.4 GHz in the B 
configuration in 2001 March, with a similar ob- 
serving strategy, except that the total observation 
time of four hours was spread out over an 8-hr 
period. 

There are many additional observations of 
Arp 299 in the VLA archives, mostly from su- 
pernova monitoring campaigns. We selected a 



^Thc Very Large Array (VLA), the Very Long Baseline Ar- 
ray (VLB A), and the Robert C. Byrd Green Bank Telescope 
(GBT) are all operated by the National Radio Astronomy 
Observatory (NRAO), which is a facility of the National 
Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 



combination of data sets covering a wide range of 
baselines to complement the new 4.9 and 8.4 GHz 
observations and to maximize the aperture-plane 
coverage of the resulting combined data. We used 
only observations with more than 200 seconds of 
integration on Arp 299. All of the selected ob- 
servations had a bandwidth of 50 MHz in each 
of two intermediate frequency channels for both 
right and left circular polarization. Each observa- 
tion used 3C 286 as the flux calibrator and was 
calibrated to the primary flux density scale of 
Baars et al. (1977), as modified very slightly by 
VLA staff to account for source variability. We 
obtained 14 sets of archival data at 8.4 GHz and 
11 sets of archival data at 4.9 GHz. We also re- 
trieved one observation at 1.4 GHz. Observations 
exist in the VLA archive at higher frequencies 
(15 and 22 GHz); however we did not use these 
because the calibration is less certain and the 
data give inconsistent results. Details of all new 
and archival observations used in this paper are 
listed in Table 1. The noise levels for each indi- 
vidual image are listed; these images were made 
with the "robust" parameter for data weighting 
that provides the best combination of high sen- 
sitivity and high resolution (Briggs 1995; Briggs, 
Schwab, & Sramek 1998). The combination of the 
archival data with the new observations provided 
very high sensitivity over a wide range of spatial 
scales at 4.9 and 8.4 GHz. Although the image 
at 1.4 GHz is not of comparable sensitivity to the 
4.9 and 8.4 GHz images, it is adequate for com- 
parisons of emission from the strongest compact 
radio sources. 

Each individual data set was calibrated, flagged, 
imaged, self-calibrated, and re-imaged using stan- 
dard image processing techniques in the Astro- 
nomical Imaging Processing System, AIPS (van 
Moorsel, Kemball, & Greisen 1996). At 4.9 and 
8.4 GHz, the multiple data sets at each wavelength 
were combined into a single data set and another 
cycle of imaging plus self-calibration and imag- 
ing produced images with high resolution and low 
noise. The data processing included the removal 
of two background radio sources within the field 
of view. They are a companion E/SO galaxy lo- 
cated 1' to the northwest (sometimes referred to as 
IC 694 (Yamaoka ct al. 1998)), and a strong back- 
ground radio galaxy ~ 1.5' southeast of Arp 299 
(centered approximately at J2000 right ascension 



ll''28'"29^ declination 58°32'05"). 

The default 8.4 GHz combined image was of 
higher resolution than the 4.9 GHz image. Thus 
the data weighting for the 8.4 GHz data was 
adjusted so that the resulting 8.4 GHz image 
matched the 0''38 x0"31 resolution of the 4.9 GHz 
image; this adjustment increased the rms noise in 
the combined 8.4 GHz image. Figures 2 and 3 
show the combined images of Arp 299 at 4.9 and 
8.4 GHz, respectively. The 1.4 GHz image (Fig- 
ures 4) could not be adjusted to match the resolu- 
tion of the 4.9 and 8.4 GHz images, so additional 
images at 4.9 and 8.4 GHz were made to match 
the resolution of the 1.4 GHz image. These im- 
ages are not shown, but are discussed occasionally 
below and are used in some of the analysis. 

All of the images show five dominant compact 
radio sources. The four strongest radio sources 
(labeled A, Bl, C, and C') have been identified 
previously by Gehrz et al. (1983). The fifth bright- 
est compact radio source, which we have labeled 
D after the Gehrz et al. (1983) notation for the 
region where this source is located, was identified 
previously at radio wavelengths by Huang et al. 
(1990) and discussed by Zhao et al. (1997). Source 
D clearly is visible in the 1.6-/im NICMOS image 
shown by Alonso-Herrcro et al. (2000), and in the 
HST 814 nm image (Figure 1). In our new ra- 
dio images. Source B has been resolved into two 
sources, Bl and B2, for the first time at radio 
wavelengths (see Figures 2 and 3); these sources 
clearly are separated in near-infrared and optical 
images shown by Lai et al. (1999) and Alonso- 
Herrero et al. (2000). 

Sources A, Bl, C, C', and D are the focus of 
this paper. 

3. VLBA Observations and Imaging 

We observed Arp 299 using Very Long Baseline 
Interferometry (VLBI) to search for possible ac- 
tive galactic nuclei (AGNs) or radio-luminous su- 
pernovae (SNe) within Sources A, Bl, C, C', and 
D. 

Three sets of VLBI observations were obtained 
in 2002 and 2003. In all cases, dual polarizations 
were observed, and Arp 299 was phase-referenced 
(Beasley & Conway 1995) to the calibration source 
J1128-h5925. The calibrator is only 0.9° away 
from the galaxy, allowing for good atmospheric 



calibration. On 2002 April 29, we observed the 
galaxy at 2.3 GHz using the VLBA^, and including 
the 100 by 110 meter GBT^ for increased sensitiv- 
ity. The observation was made with 32 MHz total 
bandwidth and 2-bit sampling, for a total recorded 
datarateof 256 Mbit s~^. One- minute scans of the 
calibrator and 2.5-minute scans of the galaxy were 
alternated; total integration time on the galaxy 
was approximately 5 hours. A second observa- 
tion was made on 2003 February 9 using only the 
VLBA. The galaxy was observed at 8.4 GHz in the 
same observing mode as in 2002; 1-minute scans 
on the calibrator and 3-minute scans of the galaxy 
were alternated. The total time on the galaxy was 
approximately 7 hours. Finally, a third follow-up 
("target of opportunity") observation was made 
on 2003 April 30/May 1 (hereafter 2003 May 1), 
using only the VLBA, at 2.3, 5.0, 8.4, and 15 GHz 
These observations used only 16 MHz total band- 
width and half the data rate of the previous ob- 
servations. Galaxy scan times ranged from 1 to 
3 minutes in the phase-referencing. The integra- 
tion time at each of the four wavelengths was se- 
lected to achieve equal noise levels at all bands; in 
the final datasets, the noises at 2.3 and 15 GHz 
were slightly higher due to ambiguous phase con- 
nection (15 GHz) and radio-frequency interference 
(2.3 GHz). Details of the VLBA observations are 
listed in Table 2. Here, the listed noise levels are 
those for which "natural" weighting was used; this 
maximizes sensitivity, at some expense in resolu- 
tion (Briggs et al., 1998). 

Calibration, imaging and self-calibration of the 
VLBA data also were carried out using AIPS. Am- 
plitude calibration was based on system tempera- 
ture measurements made during the observations, 
and has an accuracy of ^ 5%. Instrumental de- 
lays, delay rates, and phases were calibrated by 
using the phase-reference calibrator and then ap- 
plied to the galaxy. Ionospheric corrections also 
were made using data provided by the Jet Propul- 
sion Laboratory"^. The data were processed us- 
ing short correlator integration times (0.5 seconds) 
and enough spectral channels (32) to enable imag- 
ing at all locations in the merger without losses in 
peak intensity due to bandwidth or time-average 
smearing (Bridle & Schwab 1998). Self-calibration 



^lONEX files can be obtained by FTP from cd- 
disa.gsfc.nasa.gov. 



was carried out when possible in our low sig- 
nal/noise data, and was most successful for the 
first epoch 2.3 GHz observations, due to the highly 
sensitive baselines to the GBT. 

The 2.3 GHz image of Source A from 2002 (Fig- 
ure 5) contains four unresolved compact sources 
located within 0'.'15 of the VLA position for this 
source. A fifth unresolved object is detected at the 
position of Source D. The other relatively strong 
VLA sources, Bl, C, and C', were not detected. 
The initial 8.4 GHz image from 2003 February 
showed all the sources present in 2002 - four in 
Source A and one in Source D. Efforts to self- 
calibrate the data at quite low signal/noise caused 
some of these sources to disappear, while others 
appeared in the images. Therefore, we use the 
original phase-referenced images in this paper. As 
in 2002, Bl, C, and C' remained undetected. How- 
ever, a surprising result, shown in Figure 6, was 
the appearance of a new ^ 3 mJy source, dubbed 
AO (as the strongest source), only 13 mas (2.6 pc) 
from Al. This source is discussed in more detail 
in Section 5.2. 

Since we could not use self-calibration on Arp 
299 in 2003 February, we employed the check 
source for the phase referencing, Jl 127+5650 (lo- 
cated roughly 2.5° from J1128-I-5925), to estimate 
the coherence loss due to imperfect phase calibra- 
tion. Comparing the original and self-calibrated 
images for this check source, we find a coherence 
loss of 24%. The phase connection for the refer- 
ence source is good, with apparent changes (de- 
pending on baseline length) of 20° to 50° between 
successive observations of the reference source; lin- 
ear interpolation between these scans probably 
gives a phase error with respect to time of no more 
than ^ 10°, which could account for only a small 
fraction of a 24% coherence loss. Therefore, we 
infer that most of the coherence loss is caused by 
spatial structures in the atmosphere, in which case 
the phase error is likely to scale with the angular 
separation of reference source and target. Since 
Arp 299 is only 0.86° from J1128+2925, we infer 
that the loss of coherence for the galaxy should 
have been ^ 1/3 that for the check source, or 
about 8%. Therefore, we have increased the mea- 
sured source flux densities at this epoch by 8%; 
since the uncertainty in this estimate of coherence 
loss may be as much as a factor of two, we also 
have added 8% (in quadrature) to the other error 



sources in computing the final flux density errors. 

The appearance of AO motivated a target-of- 
opportunity proposal for multi-frequency VLBA 
imaging, in order to determine the spectrum of 
this new source. These observations in 2003 May 
were carried out with the VLBA alone, used lower 
bandwidth than the previous observations, and 
did not include the GBT; therefore they have 
higher rms noises than the earlier observations. 
Source AO was not detected above the 2cr level at 
2.3 GHz in May 2003. Components Al, A2, A3, 
and A4 are detected at levels consistent with their 
previous sizes and flux densities, given the much 
higher error levels. Source D changed significantly 
between 2002 April and 2003 May: the 2.3 GHz 
flux density in 2003 was more than 6 times higher 
than it was a year earlier, and Source D also be- 
came resolved at both 8 and 15 GHz (Figure 7). 

The properties and characteristics of these com- 
pact sources are discussed below. 

4. Properties of the Compact Radio Sources 
4.1. VLA Results 

In both the high-resolution images and those 
matched to the 1.4 GHz resolution, single Gaus- 
sian models were fitted to the five strongest com- 
pact radio sources. If the model fitting indicated 
that the compact source was consistent with being 
unresolved, the source was re-fitted using a point- 
source model. Results of these fits are given in 
Tables 3 and 4. Flux density errors combine the 
noise, the fitting errors, and possible systematic 
calibration offsets, added in quadrature. For well- 
calibrated data, the overall amplitude scaling rela- 
tive to Baars et al. (1977) is reported by the VLA 
Calibrator Manual^ to be 1% to 2% at 1.4, 4.9, 
and 8.4 GHz. However, because most of our data 
are archival, with a variety of different weather 
qualities and observing strategies, we have conser- 
vatively assumed flux-scale errors of 5%. 

Peak brightness temperatures at 8.4 GHz were 
calculated for Sources A, Bl, C, C, and D using 
the flux densities and sizes found at the highest 
resolution (see Table 3), and are given in Table 5. 
The brightness temperatures Tb were calculated 
using Equation (1) from Falcke et al. (1996), mod- 



*The latest version of the VLA Calibrator Manual is acces- 
sible at http://www.vla.nrao.edu/astro/. 



ified for an elliptical gaussian: 

\mJy/ \GHz/ yarcsec^/ 

(1) 
Here, S^ is the flux density at frequency i^, with 

9i and 62 being the fitted full widths at half max- 
imum of the major and minor axes of the sources. 
Total monochromatic radio powers for each source 
were derived from the matched-resolution images 
(see Table 4). These powers, given in Table 5, typ- 
ically are thousands of times the power emitted 
by the Galactic supernova remnant Cassiopeia A, 
which is assumed to have a present 8. 4- GHz fiux 
density of 460 Jy (derived from Baars et al. 1977) 
at a distance of 3.4 kpc (Reed et al. 1995). The 
inferred brightness temperatures and radio pow- 
ers are as expected from complexes of hundreds 
to thousands of supernova remnants, although we 
note that a few individual supernova remnants in 
dense environments have been observed to show 
comparable luminosities (Weiler et al. 1990; Van 
Dyk et al. 1993). 

Figure 8 shows the radio spectra of the five 
compact components listed in Table 4. On VLA 
scales, Sources A, Bl, C, and C all have spectra 
that decrease slightly with increasing frequency, 
while VLA source D has a spectrum that appears 
to peak at frequencies above 8 GHz. Two-point 
spectral indices, a {Si, oc v'^"), where 5*^ is the 
flux density at frequency v), were calculated be- 
tween 4.9 and 8.4 GHz from the images made with 
the same l'.'21x0'.'88 beam, and are listed in Ta- 
ble 5 (Supernova rates in the table are derived in 
section 5.3, below). In spite of the small frequency 
range spanned, errors in these spectral indices are 
fairly small, ^ 0.13, because the five sources under 
study are detected with high signal/noise ratios. 
The spectral indices of —0.5 to —0.65 for Sources 
A, Bl, C, and C' are consistent with collections 
of typical galactic supernova remnants (Weiler et 
al. 2002). The spectrum of Source D is more enig- 
matic; we return to it in Section 5.4. 

4.2. VLBA Results 

In all VLBA images, compact sources were vi- 
sually identified and then fitted to a gaussian using 
the AIPS task JMFIT. Details of the fits and flux 
measurements of the VLBA-scale compact sources 
are listed in Table 6 and Table 8. Random er- 



rors are assigned as the quadratic sum of the im- 
age noise levels (Table 2) and the amplitude cali- 
bration errors, estimated to be 5% at most wave- 
lengths and (conservatively) as 10% at 15 GHz. In 
addition, for data that could not be self-calibrated 
in 2003 February, the flux density at 8.4 GHz has 
been increased by 8%, with an 8% error added in 
quadrature, as inferred from observations of the 
check source J1127+5650 (see Section 3). 

An important result of this work is the VLBI 
identification of a number of compact radio 
sources with high brightness temperatures, em- 
bedded within Source A and in Source D. Ta- 
ble 7 gives the brightness temperatures of these 
milliarcsecond-scale components, measured in the 
highest sensitivity 8.4 GHz VLBI images from 
2003 February. For the unresolved sources, size 
upper limits for each dimension are taken conser- 
vatively to be half the diameter of the half-power 
beamwidth. In addition to the brightness temper- 
atures, powers relative to the galactic supernova 
remnant Cassiopeia A (Gas A) are listed, as are 
monochromatic radio powers at 8.4 GHz. Source 
D is assumed to be at the distance of Arp 299 (see 
Section 5.4); two values are given for this object, 
one when it was unresolved in 2003 February, and 
one when it was resolved in 2003 May. 

The VLBI-scale sources have very different 
spectra, as can be seen in Table 8 and in Figure 9, 
which show the individual source flux measure- 
ments at four frequencies in 2003 May. Source 
AO has an extremely inverted spectrum: in May 
2003 the 5.0 and 8.4 GHz flux densities implied 
a spectral index of a > +1.76 ± 0.37. Gompari- 
son between the (non-contemporaneous) flux mea- 
surements of a 3(7 upper limit of 0.11 mJy at 
2.3 GHz in 2002 April, and the 8.4 GHz value of 
3.23 ±0.28 mJy in 2003 February, yields a > +2.6 
for Source AO, which declines more steeply at 
long wavelengths than is expected even for syn- 
chrotron self-absorption. Source Al has a non- 
thermal spectrum, with a — —1.25 + 0.28 between 
2.3 and 8.4 GHz. Sources A2, A3, and A4 ap- 
pear to have fairly flat spectra although the S/N 
in these (short) measurements is not impressive. 
However, if we assume no variability at 2.3 and 8.4 
GHz between 2002 April and 2003 February, all 
three sources have flat or slightly inverted spectra. 
In May 2003, Source D had a relatively flat spec- 
trum, rising somewhat at the higher frequencies. 



Although we do not know definitively what spec- 
trum VLBI Source D had earlier, we do know that 
its 2.3 GHz flux density increased from 0.30 mJy 
in Feb 2002 to 2.04 mJy in May 2003 and that 
its 8.4 GHz flux density increased from 2.21 mJy 
to 3.36 mJy between Feb and May 2003, as it 
appeared to increase dramatically in size. 

5. Discussion 

5.1. Does Arp 299 host one or more 

AGNs? 

Any AGNs in Sources A, Bl, G, G' or D are 
buried in many magnitudes of extinction at visible 
and IR wavelengths (Alonso-Herrero et al. 2000). 
Sources A and Bl are thought to be the nuclei of 
the two galaxies which are merging, and both are 
associated with hard X-ray point sources (Ballo 
et al. 2003; Zezas et al. 2003), so they seem to 
be hkcly AGN candidates. [Ballo et al. (2003) 
finds hard X-ray emission and Zezas et al. (2003) 
show that it is localized at A and BL] VLBI ob- 
servations on baselines of thousands of kilometers 
are able to detect even optically obscured AGNs, 
because VLBI isolates only high brightness tem- 
perature (Tfc > 10^-10^ K) objects. In contrast, 
H II regions must have Tb < 10"*, and supernova 
remnants generally have low brightness tempera- 
tures as well. However, high brightness tempera- 
ture is not an exclusive discriminator for AGNs, 
since some young radio SNc also are detected on 
milliarcsecond scales (e.g., Bietenholz, Bartel, & 
Rupen 2001; Perez-Torres et al. 2002; Bartel & Bi- 
etenholz 2003), implying brightness temperatures 
similar to some AGNs. 

There have been a number of suggestions in the 
literature, based on observations in a variety of 
frequency bands, that source A may harbor an 
AGN (Lonsdale et al. 1992; Lai et al. 1999; Alonso- 
Herrero et al. 2000; Zezas et al. 2003). In fact, 
the total VLA-scale radio power of Source A is 
considerably higher than the cores of most low- 
luminosity AGNs (LLAGNs) and classical Seyfert 
galaxies (Ho & Ulvestad 2001; Schmitt et al. 2001; 
Thean et al. 2001; Nagar et al. 2002). Even if 
an AGN is present, it is clear from the optical 
and infrared imaging that Source A contains a 
strong star-formation complex. The inferred su- 
pernova rate in this complex, derived from radio 
flux densities and the strength of near-infrared [Fe 



II] emission, is 0.6-0.7 yr~^ (Alonso-Herrero et al. 
2000), while the VLA radio flux density imphes a 
supernova rate as high as several per year (see Sec- 
tion 5.3). Strong CO emission (Sargent & Scoville 
1991; Casoli et al. 1999) and the lack of a strong 
optical nucleus (Alonso-Herrero et al. 2000) at the 
location of source A confirm that it is heavily en- 
shrouded by dust. 

Previous VLBI measurements by Lonsdale et 
al. (1992) detected Source A with a correlated flux 
density of 1.4 mJy at 1.6 GHz on a transatlantic 
VLBI baseline, implying a source size < 1 pc with 
a radio power of 3x10^° W Hz^^ and a bright- 
ness temperature > 5x10^ K. (Smith, Lonsdale, 
& Lonsdale 1998) detected source A with corre- 
lated flux densities of 15 and 7 mJy at 1.6 GHz, on 
respective projected baselines of about 50 km and 
100 km, implying that there is considerable struc- 
ture on scales between ~ 0''005 and ^ 0''5. There- 
fore, any AGN that might be present is embedded 
in a larger region of star formation and probably 
does not make a significant contribution to the ra- 
dio power of Source A. Other star-formation com- 
plexes on similar scales have been detected with 
VLBI observations: Arp 220 (Smith et al. 1998), 
which contains an AGN but not an energetically 
significant one, and Mrk 231 (Taylor et al. 1999), 
which contains a strong AGN. 

Our detection of five compact radio sources 
with Tb > lO^K in Source A makes it unlikely that 
all the milliarcsecond-scale radio emission is AGN- 
powered, since it is highly unlikely that there are 
five AGNs within Source A. However, there could 
be a single AGN present, as inferred from the hard 
X-rays (Ballo et al. 2003; Zezas et al. 2003). In 
this regard, it is notable that source Al is the 
only milliarcsecond (<1 pc) compact object with 
an optically thin synchrotron spectrum; if there is 
an AGN present, it mostly likely is Source Al. If 
Source Al were an AGN, it would be intriguing 
to have an apparent young supernova (AO) only 
2.6 pc away from the AGN. 

The X-ray point source associated with Source 
A has a 0.1-10 keV luminosity of 4 x 10*^^ ergs s~^ 
(Zezas et al. 2003). If the X-ray emission from 
Source A were assumed to be associated with 
a single one of its milliarcsecond radio sources, 
such as Source Al, we could estimate the ra- 
dio to X-ray ratio for that source. For Source 
Al, Rx = vP^ (5 GHz)/L(2 - 10 keV) would be 



^ 2 X 10^^, consistent with values found for several 
other LLAGNs imaged on milliarcsecond scales 
(Ulvestad & Ho 2001b) and also consistent with 
the values found for the nucleui in the merger sys- 
tem NGC 3256 (Neff et al. 2003). Source Bl has 
hard X-ray emission suggestive of an AGN (Ballo 
et al. 2003; Zezas et al. 2003), somewhat weaker 
than A, but no compact AGN is detected at ra- 
dio wavelengths. Still, the radio to X-ray ratio for 
the milliarcsecond source in Bl is consistent with 
the range of values found for LLAGNs. Sources C, 
C', and D have much softer X-ray spectra (Zezas 
et al. 2003), as well as no high brightness radio 
emission, suggesting that they harbor no AGNs. 
Although the radio observations alone do not de- 
termine whether there are AGNs in Arp 299, we 
are able to say that any AGN present are unlikely 
to play an important role in the merger's overall 
energy balance. 

5.2. Young Supernovae Within Source A - 
A Supernova Factory 

Source A is particularly strong; its VLA 
8.4 GHz radio power of 1.8x10^^ W Hz^^ is more 
than 50 times the strongest non-thermal source in 
NGC 4038/9 (Neff & Ulvestad 2000) and ~27,500 
times the power of Cas A at the same wavelength. 
Source a contains at least three distinct near- 
infrared sources, possibly young star clusters, in 
its central 0.86 arcseconds (Alonso-Herrero et al. 
2000). The present VLBI imaging reveals at least 
five compact radio emitters within Source A. Four 
of these five have flat or inverted radio spectra 
between 2.3 and 8.4 GHz, a property that is char- 
acteristic of young Type II supernovae near their 
peaks (Weiler et al. 2002). Although such spectra 
also may be indicative of thermal gas associated 
with young star clusters, such clusters would have 
brightness temperatures < 10^ K, and would not 
be detected by VLBI observations. 

Source AO has an extremely inverted spectrum 
(with a > -1-1.8, Figure 9). Such an inverted spec- 
trum is characteristic of very young supernovae in 
the early days after their explosion, when the su- 
pernova blast has not yet completely penetrated 
the circumstellar mass-loss shell, which remains 
optically thick to lower frequency radio emission 
for months in "normal" SN environments, and for 
years in very dense surroundings (Chevalier 1982; 
Weiler et al. 2002). The spectra of Sources A2, A3, 



and A4 are typical of young SNe and Source Al 
has a spectrum suggesting an older SN remnant 
which has become optically thin down to frequen- 
cies of 1.4 GHz. 

If we assume that all five VLBI sources within 
A are young supernovae, they are spread over a 
total extent of 90 pc (440 mas) . The three sources 
AO, Al, and A2 are much more tightly clustered, 
with projected separations of only 2.6 pc (13 mas) 
between AO and Al, and 10.4 pc (52 mas) be- 
tween AO and A2. Assuming that these sources 
all are young supernovae, it is quite possible that 
AO and Al, and perhaps even A2, arc found within 
the same cluster of young stars. Super star clus- 
ters seen optically in other merging galaxies have 
canonical radii of ~ 4 pc, but may be as large 
as tens of parsecs (e.g., Whitmore et al. 1999). 
Alonso-Herrero et al. (2000) infer that Source A 
had the peak of its starburst 6-8 Myr ago. This 
corresponds to the main-sequence lifetime of O 
stars of ~25-30Mo (Schaller et al. 1992; Schaerer 
et al. 1993). Massive stars formed in that burst 
should just now be ending their lives on the main 
sequence, and the radio sources that we see now 
probably result from the interactions of these ex- 
plosions with their surroundings. We should there- 
fore expect to see frequent supernova explosions 
from this "supernova factory. " 

The five identified VLBI sources account for 
only 6% of the total VLA 8.4 GHz flux density of 
87 mjy in Source A, which has an optically thin 
synchrotron spectrum on a scale of hundreds of 
milliarcseconds (cf. Table 3). It is likely that most 
of the radio emission not detected on milliarc- 
second scales results from earlier supernovae that 
have now expanded and faded beyond the detec- 
tion range of VLBI and merged into the smoother 
interstellar medium. The expansion to size scales 
larger than ~ 0.2 pc typically takes a hundred 
years or so in a "normal" interstellar medium, 
while the general merging into that interstellar 
medium may take ~ 10^ yr (Condon 1992). How- 
ever, these time scales are quite uncertain in the 
denser environment of Arp 299. 

If Source A does contain five young supernovae, 
why are none detected in sources Bl, C, or G'? 
Alonso-Herrero et al. (2000) use infrared fluxes 
to infer supernova rates 0.6-0.7 yr~^ for Source 
A, but flnd supernova rates for Bl, G, and C' 
that are factors of 5-20 below that for Source A. 



We derive similar rates based on the VLA-scale 
radio emission (see Section 5.3). It seems rea- 
sonable that the number of detectable radio su- 
pernovae in B, C, and G' should therefore be no 
more than 5%-20% of the number in Source A, or 
fewer than one per region for our most sensitive 
VLBI observations. If the supernovae in A tend 
to be stronger radio emitters due to a denser inter- 
stellar environment, as we expect, the number of 
supernovae detectable in the other sources would 
be even smaller. Sources Bl, G, and G' are also 
thought to be younger starbursts (Alonso-Herrero 
et al. 2000), so they may be just beginning signif- 
icant supernova activity. 

Another alternative for the compact radio 
sources in source A is that they are X-ray bi- 
naries or microquasars. However, X-ray bina- 
ries typically have radio to X-ray ratios Rx ^ 
10"'"^ (Fender & Kuulkers 2001). Therefore, we 
would expect > 10''^ ergs s^^ in X-rays from 
Source A if the compact radio sources were X- 
ray binaries, compared to a measured value of 
~ 1.3 X 10''° ergs s"^ (Zezas et al. 2003). This 
appears to rule out active binaries as the origin of 
the milliarcsecond radio sources. 

Future long-term VLBI monitoring of Arp 299 
should enable us to test the idea that Source A 
contains a cluster of young SNe. First, new super- 
nova explosions could be discovered as they occur, 
providing a direct measure of the radio supernova 
rate and indicating whether or not all supernovae 
in Source A are likely to be radio loud. If the de- 
tected radio supernova rate were one per decade 
rather than one per year or so, much less than the 
total estimate of ^ 0.6 yr^^ (Alonso-Herrero et al. 
2000) , then this result would indicate strongly that 
only a fraction of the Arp 299 supernovae are radio 
loud, at least above the VLBA detection thresh- 
old. Second, the fading of the young supernovae 
could be measured, for comparison with observa- 
tions of young supernovae in less dense environ- 
ments; this would provide useful constraints on the 
properties of the star-forming regions, as well as 
enabling an estimate of the age of the young radio 
supernovae already detected in Source A. Third, 
we should be able to verify the nature of Source 
Al: it should fade monotonically or rise and fall 
coherently if it is a supernova, while more stochas- 
tic flux variations would suggest that it is an AGN. 



5.3. The Supernova Rate for VLA Sources 
A, Bl, C and C 

Radio sources with "normal" nonthermal spec- 
tral slopes (i.e., a < 0.0) are likely to be com- 
plexes of older supernova remnants, with the radio 
emission dominated by optically thin synchrotron 
emission. Inspection of Tables 4 and 5, where the 
resolutions at all wavelengths have been matched 
to the 1.4 GHz beam, shows that sources A, Bl, C, 
and C all have such spectra on arcsecond scales. 
The spectral indices of —0.5 to —0.65 are con- 
sistent with typical galactic supernova remnants 
(Weiler et al. 2002), allowing us to use the ob- 
served flux and the spectral index to derive the 
supernova rate in each of these arcsecond-scale 
(~100-200 pc) sources. 

There has been considerable discussion about 
the supernova rate corresponding to a given ra- 
dio flux density, with estimates varying by a fac- 
tor of 10 or so. In studying NGC 4038/9, Neff 
& Ulvestad (2000) pointed out that very differ- 
ent supernova rates will be derived, depending on 
whether the global radio emission of the merger, or 
only the emission identifiable with compact radio 
sources, was considered. This is because the bulk 
of the total radio emission of a starburst galaxy 
is caused by synchrotron-emitting electrons that 
have escaped the vicinity of the supernova rem- 
nants that accelerated them, and diffused through- 
out the galaxy. In NGC 4038 Neff & Ulvestad 
(2000) obtained the same supernova rates using 
the approach in Ulvestad (1982) and using the 
technique described by Condon (1992), as long as 
the latter approach was applied to the radio emis- 
sion from the entire galaxy. In Arp 299, as in 
NGC 4038/9 (Neff & Ulvestad 2000), we consider 
only radio emission that is compact (arcsecond 
scale or smaller), and therefore derive the super- 
nova rate for the VLA-scale sources using Equa- 
tion (7) from Ulvestad (1982) (equation 2, below). 

For Sources A, Bl, C and C, we derive the su- 
pernova rate by assuming supernova blast wave 
energies of 10^° ergs and ambient densities of 
1 cm^^. In this case, the supernova rate i?sN is 
given by: 



i?SN == 0.37/ X 



S. 



4.9 GHz 



5.0 mJyy V4885 



408 



a+0.75 



yr 



(2) 



Here, / is the fraction of the radio flux den- 
sity that can still be associated with individual 
supernova remnants on parsec scales. For the 
4.9 GHz flux densities, we use the value measured 
in the matched angular resolution images (Ta- 
ble 4). For the spectral index a between 408 MHz 
and 4.9 GHz, we use the measured value between 
4.9 GHz and 8.4 GHz (Table 5), thus assuming 
that an ensemble of supernova remnants unre- 
solved by the VLA has a net spectrum that is 
straight between 408 MHz and 8.4 GHz. The re- 
sulting estimates of the supernova rates required 
to produce the observed radio flux densities for 
sources A, B, C, and C, are given in Table 5. 

The total radio emission from source A, which 
corresponds to a power 2.7 x 10'* times that of Cas 
A, requires a supernova rate of ~ 4.7/ yr~^. The 
nonthermal sources in Bl, C, and C are consider- 
ably less powerful, with their supernova rates esti- 
mated at 0.3/-1.0/ yr^^. As stated in Section 5.2, 
the five identified VLBI sources account for only 
6% of the total flux density at 8.4 GHz, so a lower 
limit for / is ~ 0.06 for Source A. Since there may 
be other emission on scales of several parsecs that 
is too big to be seen in the milliarcsecond VLBI 
imaging, a reasonable value for / may be 0.1-0.2, 
which then would provide a calibration similar to 
those advocated by Condon (1992) and Huang et 
al. (1994). This would imply a supernova rate of 
0.5-1.0 yr^* for Source A, in good agreement with 
the result found by Alonso-Herrero et al. (2000). 

5.4. What is Source D? 

Radio source D does not appear to be in the 
main interaction region of Arp 299 (see Fig. 1 or 
2). Therefore, it is worth considering the possibil- 
ity that Source D actually is not associated with 
the galaxy, but is a background source instead, 
as suggested previously (Zhao et al. 1997). Using 
the VLA in 1980, Condon et al. (1982) reported a 
continuum radio source in the vicinity of source 
D at 1.4 GHz, with a flux density of ^4 mJy 
in a 2" beam. This flux density is much higher 
than the 1 mJy upper limit in 1989 (Huang et al. 
1990) the 0.8 mJy flux density that we derive from 
1991 archival data (Figure 5, Tables 3 and 4) the 
0.5 mJy measured in a snapshot observation on 
2003 April 29, or the VLBI upper limit from this 
work. Source D might be variable at 1.4 GHz, but 
we consider it more likely that the measurement 



by Condon et al. (1982) may have included more 
diffuse emission surrounding the source - we note 
that the VLA was in an early, hybrid, array at the 
time of these observations. 

It is possible to use the many individual ob- 
servations to investigate the possible variability 
of Source D; comparisons to the other compact 
sources can be used to infer the reality of any ap- 
parent variability. Therefore, we have measured 
the flux densities of sources A, Bl, C, C and D 
in independent data sets observed between 1989 
and 2003; results are shown in Table 9 and in 
Fig. 10 ^. We conservatively use only data from 
the high-resolution A and B configurations at 4.9 
and 8.4 GHz, where we expect diffuse radio emis- 
sion to be resolved out; this approach insures that 
measurements of the relatively isolated source D 
are not contaminated by more diffuse emission. 
Images were made with matched restoring beams 
for each array/frequency combination used, to as- 
sure that measurements were made with the same 
resolution. Sources A, Bl, C, and C were all 
found to be extended, and were fitted with single 
Gaussian sources which were consistent in size and 
orientation in all observations made with a partic- 
ular array. Source D was found to be unresolved in 
all individual observations, and was measured as a 
point source. Flux density errors were determined 
as discussed in section 4.1, and are approximately 
the size of the points shown in Figure 10. 

In a fixed configuration, sources A and Bl do 
not vary significantly at 4.9 GHz and 8.4 GHz dur- 
ing the 13 years of observations. Sources C and C' 
appear to vary, but they are very well resolved; we 
attribute their apparent variation to differing res- 
olutions and aperture plane coverages in (mostly) 
short observations. Source D clearly has changed 
in flux density at 4.9 and 8.4 GHz (Fig. 10); the 
measured flux densities in 1990 and 1991 are con- 
sistent with one another, but the more recent mea- 
surements show significantly higher flux densities. 
This trend is completey different from Sources A 
and Bl, leading us to infer that the variability of 



^Huang et al. (1990) reported a flux density of 1.8 mjy from 
the source at 8.44 GHz on 1990 March 1; this actually is a 
composite of two data sets a few days before and after that 
date, which we have reprocessed individually. Zhao et al. 
(1997) detected the source at 8.2 GHz in March 1993 (B 
array) with a flux density of 2.2±0.1 mJy; this data is in 
good agreement with nearby data we did use. 



Source D is real and is not a data-reduction arti- 
fact. 

On VLBI scales. Source D underwent a dra- 
matic flux-density increase at 2.3 GHz between 

2002 and 2003, as well as becoming significantly 
resolved and increasing in fiux at 8.4 GHz in the 
three-month interval between 2003 February 9 and 

2003 May 1 (see Figure 7). One can imagine 
several possible explanations for this change. If 
Source D is located within Arp 299, it may well 
have had a new supernova explosion in early 2003, 
or a supernova explosion that suddenly has broken 
through the circumstellar shell to become visible 
at the lowest VLBI frequency. If it is a back- 
ground quasar, it may have ejected a new "jet" 
component in a radio flare. However, one might 
expect such a flare to give rise to a much stronger 
variation at the higher frequencies, rather than in- 
creasing in flux by more than 600% at 2.3 GHz. 
Further, the apparent source size has increased 
by ~ 1 mas in three months. Even at a mod- 
est assumed redshift of z ^ 0.1, this would corre- 
spond to a large (but not unprecedented) expan- 
sion speed of fapp/c ^ 26. 

Our tentative conclusion is that Source D is 
located within Arp 299, probably within a very 
compact young star cluster. Supernova explosions 
in Source D would then account for the appar- 
ent flux-density variability seen in the multi-epoch 
VLA data, as well as for the rapid "expansion" 
and 2.3-GHz flux density increase in 2003. In sup- 
port of this conclusiton, we note that in 2001 July, 
Source D had a soft X-ray spectrum consistent 
with being a radio supernova (Zezas et al. 2003). 
This conclusion can be validated by future VLBI 
imaging of Source D. If the source is a relativistic 
jet, its size should continue to expand, whereas a 
newly emerging supernova will stay in the same 
apparent location, yielding a relatively constant 
size as a function of time. 

6. Summary 

In this paper, we have discussed the five bright- 
est radio sources in Arp 299. Four of the five 
brightest sources. A, Bl, C, and C', have VLA 
spectra characteristic of complexes of supernova 
remnants, with inferred supernova rates ranging 
from a few tenths to a few supernovae per year. 
VLBI imaging reveals five compact radio sources 
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within A, most of which have flat or inverted ra- 
dio spectra. Based on the sizes and spectra, we 
identify these as probable young supernovae or su- 
pernova remnants in a "supernova factory," with 
individual powers of hundreds to a thousand times 
that of Cas A. One compact source (AO) has an 
extremely inverted radio spectrum, with a « +2 
or greater between 2.3 and 8.4 GHz. Three of 
the compact VLBI sources are located within a 
region with a projected size of 10 pc and may well 
be within the same super star cluster or cluster 
complex. An interesting alternative possibility is 
that the one steep-spectrum source (Al) is a low- 
luminosity AGN, occurring at a separation of only 
2.6 pc (projected) from a very young supernova 
(AO). The VLBI sources include less than 10% of 
the total VLA flux density, implying that much of 
the radio power in Source A comes from older su- 
pernova remnants that have expanded and merged 
into the surrounding ISM. We do not know which 
(if any) of the separate VLBI sources in Source 
A produce the observed X-ray emission. If we at- 
tribute the bulk of the X-ray emission to one or 
more of the VLBI sources within Source A, the 
radio/X-ray ratio is far too high for all the com- 
pact objects to be beamed or unbeamed "normal" 
X-ray binaries or microquasars. Instead, the total 
radio/X-ray ratio is consistent with a mix of young 
supernovae, supernova remnants, and perhaps an 
AGN. 

Radio source D, located at a projected distance 
of 2 kpc from source A, is quite compact on VLBI 
scales, with the VLBI source apparently contain- 
ing all of the flux detected in the highest resolution 
VLA images. Its VLA 8.4 GHz and 4.9 GHz flux 
densities have both varied by factors of ^ 2 over 
the last 13 years. On VLBI scales, the Source 
D flux density at 2.3 GHz increased by a factor 
of ~ 6 between 2002 April and 2003 May, and 
its 8.4 GHz flux density increased by a factor of 
1.6 in three months in early 2003; during early 
2003, Source D also increased in size by at least 
a factor of 2 at 8.4 GHz. Although Source D 
could be a background quasar undergoing superlu- 
minal expansion, the observed properties flt more 
closely to a star-formation complex. VLBI moni- 
toring should allow us to verify the SN hypothesis 
(or not) for Source D, and a high-resolution mid- 
infrared image of Source D would help confirm its 
identity as a young cluster in Arp 299. 
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Table 1 
VLA Observations of Arp 299 Used in Combined Images 



Observing 


Frequency 


Array 


Integration 


Phase 


Derived Flux of 


Image rms 




Date 




Configurat 


ion Time 


Calibrator 


Phase Calibrator 


at Robust=0 


a 




(GHz) 




(min) 


(J2000) 


(Jy) 


(^Jy beam^^ 










1.4 GHz (20 cm) 


observations 









1991 Aug 3 


1.46 


A 


143 


J1219+4829 


0.57 


57 










4.9 GHz (6 cm) 


observations 








2000 Oct 24 


4.86 


A 


177 


J1128+5925 


0.46 


20.1 




1994 May 16 


4 


86 


BnA 


5 


J1146+5356 


0.37 


97.9 




1993 Dec 19 


4 


86 


D 


14 


J1219+4829 


0.62 


50.1 




1993 Jun 11 


4 


86 


CnB 


19 


J1219+4829 


0.60 


61.7 




1993 May 7 


4 


86 


b 


15 


J1219+4829 


0.62 


54.3 




1993 Mar 13 


4 


86 


B 


25 


J1219+4829 


0.62'' 


46.4 




1993 Jan 28 


4 


86 


BnA 


26 


J1219+4829 


0.61 


75.3 




1992 Aug 28 


4 


86 


D 


145 


J1146+5848 


0.37 


28.3 




1992 Apr 3 


4 


86 


C 


170 


J1146+5848 


0.37 


23.4 




1991 Jul 5 


4 


86 


A 


35 


J1219+4829 


0.49 


38.2 




1990 May 9 


4 


86 


A 


25 


J1219+4829 


0.62 


55.6 




1989 Jan 14 


4 


86 


A 


7 


J1146+5848 


0.44 


73.2 










8.4 GHz (3.6 cm) 


observations 








2001 Mar 19 


8.46 


B 


171 


J1128+5925 


0.28 


17.3 




1999 Apr 7 


8 


46 


D 


31 


J1219+4829 


0.61 


44.3 




1999 Jan 21 


8 


46 


C 


32 


J1219+4829 


0.68 


33.7 




1994 May 16 


8 


44 


BnA 


15 


J1146+5356 


0.34 


29.7 




1994 Fob 18 


8 


44 


D ^ A 


51 


J1219+4829 


0.66 


26.4 




1993 Dec 19 


8 


44 


D 


13 


J1219+4829 


0.65 


29.3 




1993 Jun 11 


8 


41 


CnB 


19 


J1219+4829 


0.63 


46.7 




1993 May 7 


8 


44 


B 


15 


J1219+4829 


0.64 


31.2 




1993 Mar 13 


8 


41 


B 


25 


J1219+4829 


0.66 


37.1 




1993 Jan 28 


8 


44 


BnA 


21 


J1219+4829 


0.63 


25.5 




1991 Jul 5 


8 


44 


A 


32 


J1219+4829 


0.51 


34.2 




1991 Jun 22 


8 


42 


A 


50 


J0958+6533 


1.13 


37.5 




1990 May 9 


8 


44 


A 


42 


J1219+4829 


0.62 


37.1 




1990 Mar 4 


8 


44 


A 


15 


J1217+5835 


0.45 


50.1 




1990 Feb 24 


8 


44 


A 


15 


J1217+5835 


0.47 


30.5 





^Robust — indicates that data weighting intermediate between pure natural weighting (highest sensitivity) and 
pure uniform weighting (highest resolution) was used in producing the image (Briggs 1995; Briggs, Schwab, & Sramek 
1998). 

The flux of this phase calibrator was estimated from those derived from the other archived data sets taken near 
the time of this observation. The original observation lacked an amplitude calibrator. 



Table 2 
VLBI Observations of Arp 299 



Observing 


Frequency 


Array 


Data 


Integration 


Image rms 


Beam Size 


Date 




Configuration 


Rate 


Time 


at Robust=5'' 


at Robust— 5^ 




(GHz) 




(Mbit s-'-) 


(min) 


(A'Jy 


beam ) 


(mas) 


2002 Apr 29 


2.27 


VLBA+GBT 


256 


330 




35 


5.58 X 4.49 


2003 Feb 09 


8.41 


VLBA 


256 


438 




51 


1.84 X 1.27 


2003 May 01 


2.27 


VLBA 


128 


58 




162 


6.19 X 4.64 


2003 May 01 


4.98 


VLBA 


128 


53 




129 


2.56 X 1.87 


2003 May 01 


8.41 


VLBA 


128 


53 




137 


1.66 X 1.26 


2003 May 01 


15.36 


VLBA 


128 


167 




152 


0.83 X 0.67 



''Robust — 5 indicates that pure natural weighting of the data was used in producing the images (Briggs 
1995; Briggs, Schwab, & Sramek 1998). 
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Table 3 
High Resolution Measurements of the Five Compact VLA Sources 



Source 


RA 


Dec 


S" (1.4 GHz) 


S" (4.9 GHz) 


S" (8.4 GHz) 


Major Axis^ 


Minor Axis^ 


P.A.'' 




(J2000) 


(J2000) 


(mjy) 


(mJy) 


(mJy) 


(arcscc) 


(arcscc) 


(deg) 


A 


11 28 33.63 


58 33 46.67 


163.3±8.2 


101.1±5.0 


77.2±3.9 


0.38 


0.30 


127 


Bl 


11 28 30.99 


58 33 40.78 


41.6±2.1 


11.1±0.6 


9.1±0.5 


0.27 


0.25 


100 


C 


11 28 30.65 


58 33 49.29 


21.8±1.1 


6.9±0.4 


5.9±0.3 


0.83 


0.74 


32 


C' 


11 28 31.33 


58 33 50.02 


12.2±0.3 


4.7±0.3 


4.1±0.2 


0.49 


0.36 


28 


D'' 


11 28 33.01 


58 33 36.55 


0.91±0.08 


1.94±0.10 


2.43±0.12 


0.17 


0.11 


168 








Image 


Properties 
















1.4 GHz 


4.9 GHz 


8.4 GHz 








Restoring Beam 






l"21x0"88 


0"38x0"31 


"38x0"31 








Beam P. A. 






-38° 


17° 


17° 








rms (^Jy beam^ ) 






67 


17 


13 









^Flux densities derived from Gaussian fits to the radio components 

•^Measured at 8.4 GHz 

■^Source D is unresolved except possibly at 8.4 GHz. 

Table 4 
Matched Beamsize Flux Densities of the Five VLA Sources 



Source 5(1.4 GHz) 
(mJy) 



S'(4.9 GHz) 5(8.4 GHz) 
(mJy) (mJy) 



A 


163.3±8.2 


115.0±5.8 


87.1±4.4 


Bl 


41.6±2.1 


20.1±1.0 


14.5±0.7 


C 


21.8±1.1 


11.0±0.6 


7.7±0.4 


C' 


12.2±0.6 


6.7±0.3 


5.1±0.3 


D 


0.91±0.08 


1.73±0.09 


2.45±0.32 



Note. — All measurements made from images with 
1.21 X 0.88 beam, Position Angle -38 degrees E of N. 



Table 5 
Derived Physical Properties of the VLA Sources A, Bl, C, C, D 



Source 


"5-8 


Tb (8.4 GHz) 


P(8.4 GHz) 


P/Pca=A 


SN Rate 






(10'^ K) 


(W Hz-i) 




(y-') 


A 


-0.51 ±0.13 


17.1 


1.8 X 10^^ 


27,500 


0.5-1.0 


Bl 


-0.60 ±0.13 


3.4 


2.9 X 10^^ 


4580 


0.1-0.2 


C 


-0.65 ±0.13 


0.2 


1.6 X 10^1 


2430 


0.05-0.1 


C' 


-0.50 ±0.13 


0.6 


1.0 X 10^^ 


1610 


0.03-0.06 


D 


+0.64 ±0.13 


3.3 


4.9 X 10^° 


770 





Note. — Spectral indices, total radio powers, and supernova rates are measured 
from the matched-resolution data given in Table 4. Peak brightness temperatures 
are measured from the highest resolution data, given in Table 3. Supernova rates 
are derived using Equation (2) in Section 5.3, witn /— 0.1— 0.2. No supernova rate 
is estimated for Source D because its spectrum docs not fit the assumptions made 
in deriving Equation 2. 
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Table 6 
VLBI Component Flux Densities 



Source 


RA 


Doc 


5(2.3 GHz)" 


S(8.4 GHz)" 


Major Axis 


Minor Axis 


P.A. 


Designation 


(J2000) 


(J2000) 


(mjy) 


(mJy) 


(mas) 


(mas) 


(dog) 


AO 


11 28 33.6212 


58 33 46.707 


< 0.11 


3.23 ±0.28 


< 0.92 


< 0.64 






Al 


11 28 33.6199 


58 33 46.699 


1.91 ±0.19 


0.60 ±0.09 


< 0.92 


< 0.64 






A2 


11 28 33.6218 


58 33 46.655 


0.40 ±0.05 


0.67 ±0.09 


<0.92 


< 0.64 






A3 


11 28 33.5931 


58 33 46.560 


0.30 ±0.05 


0.39 ±0.08 


< 0.92 


< 0.64 






A4 


11 28 33.6502 


58 33 46.538 


0.35 ±0.05 


0.40 ±0.08 


<0.92 


< 0.64 






D (2002 Apr 29) 


11 28 33.0108 


58 33 36.549 


0.30 ±0.03 




< 2.79 


< 2.25 






D (2003 Fob 9) 


11 28 33.0108 


58 33 36.549 




2.21 ±0.20 


< 0.92 


< 0.64 






D (2003 May 1) 


11 28 33.0107 


58 33 36.549 


2.04 ±0.26 


3.36 ±0.36 


1.8 


1.1 


141 



"Except as noted, 2.3 GHz measurements are from 2002 April 29 and 8.4 GHz measurements are from 2003 February 9 



Table 7 
Derived Physical Properties of the VLBI Sources 



Source 




a2-s'' 


(K) 


P (8.4 GHz) 
(W Hz^i) 


P/P(CasA) 


AO 




> 2.58 ±0.26 


> 1.4 X 10** 


6.5 X 10^° 


990 


Al 




-0.88 ±0.14 


> 2.6 X 10^ 


1.2 X 10^° 


180 


A2 




+0.39 ±0.14 


> 2.9 X 10^ 


1.3 X 10^° 


200 


A3 




+0.20 ±0.20 


> 1.7 X 10^ 


7.8 X 10^^ 


120 


A4 




+0.10 + 0.19 


> 1.7 X 10^ 


8.0 X 10^^ 


120 


D (2003 Feb 


9) 




> 9.5 X 10^ 


4.3 X 10^° 


680 


D (2003 May 


1) 


+0.38 + 0.13 


4.3 X lO' 


6.7 X 10^° 


1040 



'^Spectral indices between 2.3 and 8.4 GHz for all eomponents within Source A 
are computed using the lowest noise data from 2002 April (2.2 GHz) and 2003 
February (8.4 GHz), since there is no evidence of variability in these components. 
Source D appears to have varied, so only the simultaneous data from 2003 May 1 
were used for computing its spectrum. 

Brightness temperatures for Al, A2, A3, and A4 are based on 8.4 GHz mea- 
surements from 2003 February; two values are given for Source D, for two epochs. 



Table 8 
VLBI Flux Densities at Epoch 2003 May 01 



Source 


S(2.3 GHz) 


S(5.0 GHz) 


S(8.4 GHz) 


S(15 GHz) 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


AO 


< 0.49 


1.08 ±0.16 


2.72 ±0.19 


2.61 ±0.38 


Al 


2.10 + 0.26 


0.77 ±0.15 


0.41 ±0.14 


< 0.46 


A2 


0.71 ±0.17 


0.94 ±0.16 


0.45 ±0.14 


< 0.46 


A3 


< 0.49 


0.45 ±0.14 


< 0.41 


< 0.46 


A4 


< 0.49 


0.50 ±0.14 


0.27 ±0.14 


< 0.46 


D 


2.04 ±0.26 


1.78 ±0.22 


3.36 ±0.36 


3.45 ± 0.54 



16 



Table 9 
Individual-Epoch Flux Densities of VLA Sources 



Date 


Config. 


A 


Bl 


C 


C' 


D 






(mjy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


4.9 GHz 


2003 May 01 


vlba 


3.74 ±0.34 


< 0.39 


< 0.39 


< 0.39 


1.78 ±0.16 


2002 Apr 29 


A 


101.7 ±5.1 


10.95 ± 0.59 


4.50 ±0.49 


3.69 ±0.32 


1.54 ±0.13 


2000 Oct 24 


A 


103.7 ±5.2 


11.90 ±0.60 


7.18 ±0.37 


4.56 ±0.23 


2.15 ±0.11 


1994 May 16 


BnA 


107.2 ±5.4 


14.54 ±0.73 


7.23 ±0.37 


4.05 ±0.21 


1.42 ±0.07 


1993 May 07 


B 


117.6 ±5.9 


20.26 ± 1.03 


11.02 ±0.59 


6.46 ±0.37 


1.11 ±0.10 


1993 Mar 13 


B 


113.5 ±5.7 


18.82 ±0.94 


10.74 ±0.54 


5.38 ±0.28 


1.40 ±0.08 


1993 Jan 28 


BnA 


103.8 ±5.2 


15.16 ±0.76 


7.34 ±0.38 


4.65 ±0.24 


1.29 ±0.07 


1991 Jul 05 


A 


99.5 ±5.0 


10.86 ±0.55 


4.08 ±0.24 


4.17 ±0.23 


1.17 ±0.07 


1990 May 09 


A 


99.1 ±5.0 


10.91 ±0.55 


5.15 ±0.30 


3.36 ±0.19 


1.53 ±0.08 


1989 Jan 14 


A 


99.6 ±5.0 


10.41 ±0.56 


3.29 ±0.38 


2.78 ±0.26 


1.31 ±0.11 


8.4 GHz 


2003 May 01 


VLBA 


4.05 ±0.41 


< 0.41 


< 0.41 


< 0.41 


3.36 ±0.36 


2003 Feb 09 


VLBA 


5.29 ±0.33 








2.05 ±0.12 


2002 Apr 29 


A 


69.58 ±0.35 


5.90 ±0.30 


0.66 ±0.07 


1.39 ±0.07 


2.23 ±0.11 


2001 Mar 19 


B 


82.93 ±4.15 


12.05 ±0.60 


6.56 ±0.33 


4.62 ±0.23 


2.52 ±0.13 


1999 Oct 28 


B 


78.06 ±3.90 


11.06 ±0.57 


6.04 ±0.36 


3.30 ±0.21 


2.18 ±0.13 


1994 May 16 


BnA 


78.93 ± 3.95 


7.77 ±0.39 


3.47 ±0.19 


2.73 ±0.15 


2.35 ±0.12 


1994 Feb 18 


A 


74.66 ± 3.73 


7.21 ±0.37 


1.62 ±0.15 


2.82 ±0.18 


2.54 ±0.13 


1993 May 07 


B 


84.04 ±4.20 


11.67 ±0.60 


6.54 ±0.38 


3.85 ±0.23 


1.77 ±0.11 


1993 Mar 13 


B 


80.30 ± 4.01 


11.17±0.56 


5.07 ±0.26 


3.61 ±0.19 


2.21 ±0.11 


1993 Jan 28 


BnA 


74.26 ± 3.71 


7.55 ±0.38 


3.30 ±0.19 


2.54 ±0.14 


1.91 ±0.10 


1991 Jul 05 


A 


70.27 ±3.51 


6.21 ±0.32 


1.68 ±0.17 


1.85 ±0.15 


1.35 ±0.08 


1991 Jun 22 


A 


74.05 ± 3.70 


6.69 ± 0.34 


2.41 ±0.20 


2.59 ±0.17 


1.58 ±0.09 


1990 May 09 


A 


74.15 ±3.71 


7.22 ±0.37 


1.97 ±0.14 


2.64 ±0.16 


1.71 ±0.09 


1990 Mar 04 


A 


71.35 ±3.57 


6.61 ±0.34 


1.48 ±0.15 


1.91 ±0.13 


1.50 ±0.08 


1990 Feb 24 


A 


73.83 ±3.69 


6.81 ±0.35 


1.12 ±0.13 


1.88 ±0.13 


1.57 ±0.09 
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Fig. 1.— FS'rWFPC2 814 nm image, made from 
public archival data. North is up, East is left. The 
two components of the merger are clearly distin- 
guishable. Radio source A is located in the center 
of the eastern galaxy; its IR luminosity is ~50% of 
the total IR luminosity of the galaxy pair. Sources 
B and C are in the central region of the western 
galaxy, B lower and C upper. This image shows 
approximately the same field of view as the radio 
images in Figures 2-6. 
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Arp 299 



4.9 GHz (6 cm) 




1128 34.0 33.0 32.0 31.0 

Right Ascension (J2000) 



Fig. 2.— 4.9 GHz (6 cm) VLA image of Arp 299. 
Radio sources A, B, C, C, and D are indicated. 
Contour levels are -0.065, 0.065, 0.120, 0.2, 0.4, 
0.7, 1, 2, 4, 7, 10, 20, and 40 mJy beam"!, the 
peak flux density is 46.0 mJy beam^^, the rms 
noise is 15.7 /iJy beam~^, and the restoring beam 
size is 0''38 x 0''31 with position angle of 17°. 
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Fig. 3.— 8.4 GHz (3.6 cm) VLA image of Arp 299. 
Contour levels are -0.070, 0.070, 0.121, 0.200, 
0.400, 0.700, 1,2,4, 7, 10, 20, and 39 mJy beam^^, 
the peak flux density is 39.1 mJy beam^^, and 
the rms noise is 14.1 /xJy beam^^. The data have 
been weighted to match the 4.9 GHz image shown 
in Figure 2, with a restoring beam of 0''38 x 0''31 
with position angle of 17°. 



Arp 299 



1.4 GHz (20 cm) 




1128 34.0 33.0 32.0 31.0 

Right Ascension (J2000) 



Fig. 4.— 1.4 GHz (20 cm) VLA image of Arp 299. 
Contour levels are -0.27, 0.27, 0.4, 0.7, 1, 2, 4, 7, 
10, 20, 40, 70, and 100 mJy beam"\ the peak 
flux density is 101.2 mJy beam^^, the rms noise is 
69.2 /iJy beam~^, and the restoring beam size is 
i;'21 x0"88 in position angle of -38°. Extensive 
diffuse emission is clearly present, even in these A 
configuration observations. 
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Arp 299 
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Fig. 5. — 2.3 GHz image of Source A, using 
the VLBA together with the GBT, from data 
acquired on 2002 April 29. Contour levels are 
at (—3, 3, 6, 24) times the rms noise of 3.7 x 
10~^ mjy beam~^, and the peak flux density 
is 1.76 mJy beam^^. The coordinate of (0,0) 
mas corresponds to (a, <5) = (ll''28'"33.62201", 
58°33'46'.'6100). VLBI sources Al, A2, A3, and 
A4 are indicated. The cross marks the position of 
source AO, not detected at 2.3 GHz. The beam 
size is 5.58 x 4.49 mas in position angle —28°. 





: 1 .76 m Jy/beam 
Uvs = 3.7E-05 ' (-3, 3, 6, 24) 



20 10 -10 

milliarcsec 
Peak = 3.04 mjy/beam 
Leva = 5.1 E-05 ' (-3, 3, 6, 24 ) 





Fig. 7. — VLBA images of source D from 2003 
May 1. The previously (February 2003) unre- 
solved Source D is clearly extended in these obser- 
vations, with a size of 0.35 x 0.22 pc. The left panel 
is an 8.4 GHz (3.6 cm) image, with a peak flux den- 
sity of 1.71 mJy beam~^; contour levels are at (—3, 
3, 6) times the rms noise of 0.137 mJy beam~^. 
The right panel is a 15 GHz (2 cm) image with 
a peak flux density of 1.53 mJy beam~^; contour 
levels are at (—3, 3, 6) times the rms noise of 0.152 
mJy beam~^. 



Fig. 6. — (The central region of Source A, showing 
the appearance of Source AO in the 2003 February 
observations. The left panel is a 2.3 GHz VLBA 
image from 2002 April 29, using the same contour 
levels as figure 5: (—3, 3, 6, 24) times the rms 
noise of 3.7 x 10~^ mJy beam~^ The right panel is 
an 8.4 GHz VLBA image of the same region, from 
2003 February 9. Contour levels are at (—3, 3, 6, 
24) times the rms noise of 5.1 x 10~^ mJy beam"-'^, 
and the peak flux density is 3.04 mJy beam~^. 
The 8.4 GHz image shown here has been restored 
with a beam identical to that of the 2.3 GHz im- 
age, to facilitate comparison. The cross in the 
2.3 GHz image (left panel) indicates the location 
of component AO. 
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Fig. 8. — Radio spectra of VLA sources A, Bl, C, 
C, and D. The measurements shown here are from 
the matched-resolution combined images, all with 
a restoring beam of 1''21 x0''88. Sources A, B, C, 
and C all have spectra consistent with Galactic 
supernova remnants. The spectrum of Source D 
appears to peak above 8 GHz. 




Fig. 9. — Radio spectra of the VLBI sources. 
These measurements are from the simultaneous 
observations done on 1 May 2003. Sample error 
bars at log(S/mJy)= —0.5 and 0.0 are shown on 
the right-hand side of the plot. Error bars at log 
(S/mJy)=: 0.5 are smaller than the sizes of the 
points. Source AO has a spectrum characteristic 
of a very new supernova, A2, A3, and A4 are all 
consistent with young supernovae in dense envi- 
ronments, and Al is suggestive of an older super- 
nova or supernova remnant or an AGN. 
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Fig. 10. — Flux-density history of high- resolution 
observations of three radio sources in Arp 299, as 
observed with resolution of ^ 1" or better with the 
VLA at 4.9 and 8.4 GHz. Milliarcsecond measure- 
ments with the VLB A at 5.0 and 8.4 GHz also are 
included for Source D, which appears completely 
unresolved at these frequencies. The different 
VLA configurations and the VLBA are indicated 
by different symbols, and typical flux-density er- 
rors are approximately the same size as the sym- 
bols. Note that Sources A and Bl show increased 
flux densities as the resolution is degraded, but 
appear to show no significant variability over time 
for a given VLA configuration. The unresolved 
Source D, in contrast, shows significant time vari- 
ability of its flux density, and no systematic trends 
for higher flux density at poorer resolution. 
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